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ABSTRACT 

Interest in the spatial distribution of dark matter (DM) velocities in galaxy clusters has grown 
recently in light of the improved capability of determining its degree of anisotropy, /3, from (either 
separate or joint) analysis of several different sets of measurements. Since cluster evolution is a highly 
non-linear hierarchical process, detailed theoretical expectations for /3 and its profile can only be 
obtained from numerical simulations. We report statistical results for /3 from a sample of some 6000 
cluster-size halos (at redshift zero) identified in a ACDM hydrodynamical adaptive mesh refinement 
simulation done with the Enzo code. These include profiles of j3 in clusters with different masses, 
relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. 
Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos 
(across elliptical shells extending to at least ~ 1.5r v j r ), universal patterns are found when these are 
averaged over halo masses, redshifts, and relaxation stages. These are characterized by a very small 
velocity anisotropy at the halo center, increasing outward and leveling off at ~ 0.1 — 0.2 of the virial 
radius in lower mass and redshift halos. We also find that at radii larger than about 0.2r v j r , j3 tends to 
be lower in spherical than in elliptical halos. This finding may help in sharpening the contrast between 
cluster shape measurements (e.g., by the CLASH project) and results from simulations. Our analysis 
does not indicate that there is significant correlation (found in some previous studies) between the 
radial density slope, 7, and /? at large radii, 0.3 r v - a < r < r vlI . 
Subject headings: Methods: Numerical - Galaxies: clusters: general 



1. INTRODUCTION 

Dark matter (DM), the main mass constituent of 
galaxy clusters, dominates the dynamics of intracluster 
(IC) gas and member galaxies. The DM mass density 
profile was until recently the only cluster property that 
could be inferred from simulations and tested against 
observational data. Some effort is now devoted to de- 
termine also the DM velocity anisotropy either by us- 
ing the gas temperature as a tracer of the DM velocity 
anisotropy, a method which is applicable at intermediate 
radii (Host et al. 2009), or by examining galaxy veloci- 
ties, as has recently been demonstrated in the analysis 
of A1689 measurements (Lemze et al. 2011). It is in fact 
our plan to apply the latter procedure to 14 additional 
relaxed X-ray clusters in the CLASH program (Postman 
et al. 2011). 

N-body simulations (for various cosmological models) 
suggest a nearly universal velocity anisotropy profile 
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(Cole & Lacey 1996; Carlberg et al. 1997; Colin, Klypin, 
& Kravtsov 2000; Diemand, Moore, & Stadel 2004; Ra- 
sia et al. 2004; Wojtak et al. 2005), similarly to the 
universal DM density profile deduced from simulations 
(Navarro, Frenk, & White 1997, hereafter NFW; Moore 
et al. 1998) and various observations (X-ray: e.g. Pointc- 
couteau, Arnaud, & Pratt 2005; Vikhlinin et al. 2006; 
Schmidt & Allen 2007; Arnaud, Pointecouteau, & Pratt 
2008; galaxy velocity distributions: Diaferio, Geller, & 
Rines 2005; SZ measurements: Atrio-Barandela et al. 
2008; strong and weak lensing measurements: Broad- 
hurst et al. 2005a, hereafter B05a; Broadhurst et al. 
2005b, hereafter B05b; Limousin et al. 2007; Medezinski 
et al. 2007; Lemze et al. 2008, hereafter L08; Broadhurst 
et al. 2008; Zitrin et al. 2009, 2010, 2011; Umetsu et al. 
2010). If both the density and velocity anisotropy profiles 
are indeed universal, they must be correlated. Hansen & 
Moore (2006, hereafter HM06) have recently argued for 
a universal relation between the DM radial density slope 
j(r) and the velocity anisotropy /3(r) for structures at 
virial equilibrium. Their deduced relation was claimed 
to hold for various systems, including disk galaxy merg- 
ers, simulated halos undergoing spherical collapse, and 
CDM halos both with and without cooling. 

However, while an analysis of 6 high-resolution sim- 
ulated galactic halos from the Aquarius project, carried 
out by Navarro et al. (2010), exhibited a reasonably good 
fit to the HM06 relation in the inner regions, large de- 
viations were reported outside r_2, the radius at which 
the profile slope reaches —2. Analogous results were ob- 
tained in a study conducted by Tissera et al. (2010), in 
which they resimulated 6 (Aquarius) galactic halos, con- 
structed so as to include metal- dependent cooling, star 
formation, and supernova feedback. In 3 of the halos a 
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rather good match to the HM06 relation was found at 
small radii, 2 kpc ■ hr x < r < r_2, but no corresponding 
match was found in the other 3 halos. No evidence is 
seen for the HM06 relation at large radii, r > r_2, in 
any of the six clusters. Such a relation between the DM 
density and velocity anisotropy is of interest for both fun- 
damental (HM06) and practical reasons, since the latter 
quantity is not easily measurable, whereas the density 
profile can be determined in several different ways based 
on different sets of measurements. 

We report the results of an analysis of a large num- 
ber of cluster-size halos drawn from an Adaptive Mesh 
Refinement (AMR) cosmological simulation. The large 
number of halos at different redshifts allows us to address 
the dependence of the DM velocity anisotropy profile on 
redshift, halo mass, degree of relaxation, modeled both 
as spherical and triaxial DM configurations, and to ad- 
dress also the -f-f3 relation. The outline of the paper is as 
follows. In §[5] we describe the simulation dataset, and 
in § [3] we describe how we infer the radial profiles of the 
density and velocity anisotropy in spherical and elliptical 
shells. In §|4]we specify our criteria for relaxed halos, and 
in § [5] we present the j3 profiles for different halo mass, 
redshift, and relaxation stages, and the deduced 7-/3 
relation. We conclude with a summary in § [BJ 

2. THE SIMULATION 

Clusters of galaxies were drawn from a cosmological 
AMR simulation performed with the hydrodynamical 
ENZO code developed by Bryan & Norman (1997; see 
also Norman & Bryan 1999; Norman et al. 2007), as- 
suming a spatially flat ACDM model with the parame- 
ters n m = 0.3, n b = 0.04, Qcdm = o.26, n A = 0.7, 

h = 0.7 (in units of 100 km/s/Mpc), and cr 8 = 0.9. The 
hydrodynamics in the AMR simulation used an ideal gas 
equation of state (i.e., neither radiative heating, cooling, 
star formation or feedback were included), with a box 
size of 512 h _1 Mpc comoving on a side with 512 3 DM 
particles, and DM mass resolution of about 10 11 tig 7 M Q . 
The root grid contained 512 3 grid cells, and the grid was 
refined by a factor of two, up to seven levels, providing 
a maximum possible spatial resolution of 7.8 (1 + z) 
h _1 kpc (this resolution is dependent on the criteria for 
refinement of the adaptive mesh, and we used the actual 
resolution when analyzing the halos). For more details 
on the simulation setup and analysis, see Hallman et al. 
(2007), in particular Section 2.2. Work reported here is 
based on analysis of halos found using the HOP halo- 
finding algorithm (Eisenstein & Hut 1998). 

To find the desired halo DM properties we extracted 
particle positions and velocities from the raw data. Par- 
ticles within a cube with comoving side of 16 h _1 Mpc 
were extracted. This ensured that both the halo and 
a sufficiently large surrounding region was available for 
examination. 

3. RADIAL PROFILES 

Radial profiles were extracted in both spherical and tri- 
axial shells. The mass distribution is described in terms 
of the axial ratios of the density surface contours. As- 
suming that the density distribution is stratified in sim- 
ilar ellipsoids, it is possible to determine the axial ra- 
tios without knowledge of the radial density distribution 



(Dubinski & Carlberg 1991, but see also other works e.g. 
Katz 1991; Warren et al. 1992; Jing et al. 1995; Jing & 
Suto 2002; Allgood et al. 2006, and references therein). 
The mass density in a triaxial configuration, p = p(r e ), is 
specified in terms of the elliptical distance in the eigen- 
vector coordinate system of the halo particles, r e , 

where q and s are the normalized axial ratios with s ^ 
These ratios can be derived from the tensor 

«ii=Ef (2) 
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through 
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where the sum is over all the particles, and M xx , M yy , 
and M zz are the principal components of the diagonal- 
ized tensor, with M xx ^ M yy ^ M zz . An advantage of 
this scheme is the equal weighting given to each particle 
irrespective of its radial position. The large number of 
particles in each halo allows accurate determination of 
the axial ratios. In practice, the value of r e in My is 
not known in advance, due to its dependence on q and s 
(which we want to determine) through eq. [1] The axial 
ratios are therefore determined iteratively. My is initially 
calculated assuming that the contours arc spherical, so 
that q = s = 1. Particle positions are first rotated into 
the diagonalized frame of My , where only particles inside 
of the ellipse volume were taken (a sphere, in the first it- 
eration) . The values of q and s are determined from My 
and then used to recalculate r in this new frame and fed 
back into the My relation to determine iterated values of 
q and s. When the input values match the output values 
within a certain tolerance, convergence to the true axial 
ratios is achieved. 

In each iteration new values for q and s are determined, 
so the halo volume is deformed. We kept the magnitude 
of the longest axis equal to 2r vir of the original spherical 
radius. This radius was taken since on the one hand its 
volume contains a large number of particles, > 10 3 , and 
on the other hand at these radii the ellipticity is quite 
constant, see fig. [T] Thus, during the volume deformation 
only the two smaller axes were changed. In this figure we 
plotted the averaged q and s values over all halos at z = 
at different portions of the virial radius for the main 
axis. The halo ellipticity first decreases a small amount, 
until ~ (1.5 — 2)r V j r , in agreement with Allgood et al. 
(2006) who found that halos become more spherical up to 
r = r v i r , then their ellipticity increases. In other words 
on average the halos are more elliptical at small radii 
then with increasing radius they become more spherical, 
and then elliptical again. However, the change is small 
and over the radius range (0.3 — 3)r v i r the ellipticity is 
quite constant, < q >« 0.66 and < s 0.5, especially 
compared to the large scatter. In addition, part of the 
increasing ellipticity at large radii is due to the presence 
of infalling halos. 

The DM velocity anisotropy profile for each halo was 
determined as follows: We first identified the halo center 
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parameters. The best fit was found by minimizing 
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Fig. 1. — Profiles of the average axes ratio at z = 0. Plotted 
are the average values of q (solid curve) and s (dashed curve), over 
all halos at volumes with a different main axes value at z = 0. 
The uncertainty is calculated such that Aqi = qi/\J~Nl and As^ = 
Si/^/Ni when qi and s; are the axes ratio of halo i, and Ni is 
the particle number inside the relevant ellipse. We checked for the 
most massive halo and found that indeed in the particle number 



range of 10 2 



10 the fractional error in the derived axes ratio is 



with the peak of the surrounding 3D density distribution 
and then determined the proper (non-comoving) veloci- 
ties of the DM particles with respect to the cluster cen- 
ter by subtracting the velocity of the halo center. This 
procedure was carried out for 15 equally spaced shells 
within the virial radius, a division that yields DM par- 
ticle counts of the same order of magnitude in each bin. 
Logarithmic spacing was impractical due to the low spa- 
tial resolution. The DM velocity anisotropy in each shell 
was calculated as 



= 1 



2al 



(4) 



where ay, erg, and er^ denote the radial, polar, and az- 
imuthal velocity dispersions, respectively. Shells contain- 
ing less than 10 DM particles were excluded by virtue of 
their statistical insignificance. For example, in a spher- 
ical halo and for all halos a total of 58 such shells were 
present at z = 0. We only considered halos containing 
at least 10 3 particles, so as to obtain robust results inde- 
pendent of numerical artifacts (as has also been done by 
Neto et al. 2007). Since our DM mass resolution is ap- 
proximately 10 11 h ( ^"y M , we examined all halos having 
M vir > 10 14 hoj M . 

For constructing the DM density profile we used the 
same binning and halo center definition as in DM veloc- 
ity anisotropy profile, and averaged the DM density over 
spherical shells. The radial density slope is defined as 



7M 



dlog[puM(r)} 
dlog[r] 



(5) 



For comparison with the radial density slope derived 
from a fit to an NFW profile, we fitted the resulting dis- 
tribution to an NFW profile, p. 



NFW 



(ri/r s )(l+n/r s ) 2 ' 

where r s and p s are a scale radius and the density at this 
radius, respectively, both of which were treated as free 



N bins 

X^Y, H( P >)-\og(pr W (P0,r s )f , (6) 

i=l 

where each bin was assigned equal weight. 

4. CRITERIA FOR RELAXED CLUSTERS 

The distinction between relaxed and unrelaxed clusters 
was made according to criteria laid down by Thomas et 
al. (2001) and Neto et al. (2007). These address (i) the 
displacement between the center of mass r cm and the 
potential minimum r p , and (ii) the virial ratio 2T/\U\. 
For the first criterion we defined a normalized offset, 
Soffset = \r p - r cm \/r vir (see also Duffy et al. 2008); for 
the second criterion we computed the total kinetic and 
gravitational energies of the halo particles within r V i r . 
When halos were modeled as triaxial, their major axes 
were set equal to the virial radii of the respective spher- 
ical configurations. All relevant calculations were per- 
formed for particles lying within the virial radius. For 
the estimation of T we subtracted the motion of the halo 
center, whereas U was calculated using a random sam- 
ple of 1000 particles. We controlled the precision level 
of this method by (a) repeating the calculation 10 times 
for the most massive halo (the one containing the largest 
number of DM particles) , which generated a relative dif- 
ference of (1.4 ± 1)%, and (b) calculating U in a single 
halo, using 10 4 particles. The relative average difference 
produced by this method was 0.8%. 

In equilibrium s ffset would be expected to vanish, and 
the virial ratio would approach a value slightly higher 
than unity, since even in relaxed systems there always 
is some infalling matter. While the two criteria are re- 
lated to the degree of relaxation in a straightforward 
manner, the boundary levels between the two phases are 
quite arbitrary. For example, Neto et al. (2007) adopted 
Soffsct = 0.07 and 2T/\U\ = 1.35. 

A possible third criterion to distinguish between re- 
laxed and unrelaxed systems is the substructure level, de- 
fined here as the displacement between the density peak 
r<i and the center of mass r cm , with the latter quan- 
tity calculated using particles within the virial radius. 
The displacement was normalized with respect to the 
virial radius, s su b = \rd — r cm \/r V i r . Since s su b and 

^offset 

turned up to be strongly correlated (R — 0.8), this addi- 
tional criterion was employed only once for comparison. 

5. RESULTS 

5.1. The DM velocity anisotropy profile 

Average DM velocity anisotropy profiles of all high- 
mass halos, M v j r ^ 10 14 h^ 7 Mq, are plotted in fig- 
ure^ (including la uncertainty regions) for different red- 
shifts. The profiles are essentially similar at small radii, 
r < 0.3r v i r , roughly independent of redshift. However, at 
larger radii, r ~ 0.7r V j r , values of (3 are somewhat higher 
in high redshift halos. At even larger radii, r > r vlr , /3 
is lower at higher redshifts. As the redshift increases, 
the scatter in /3 increases as a function of radius from 
~ 0.3r vlr . In figure [3] we illustrate the DM velocity 
anisotropy profiles for two mass ranges at two redshifts. 
The 100 most and least massive halos are compared at 
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Fig. 2. — Velocity anisotropy profiles at different redshifts. The 
relatively small variation of the anisotropy in the region 0—1 r v i r is 
shown separately in the upper panel. The central curves represent 
the mean value, with the the ± 1— er uncertainty around the mean 
is marked by the upper and lower curves. 

rcdshift z = 0, and a similar comparison is made at red- 
shift z = 2 for the 10 most and least massive halos. As 
is apparent from the plot, at z = the beta profile of 
the high-mass halos increases with radius; no such ten- 
dency is visible in the low-mass range. This result is 
consistent with the behavior seen in figure [5] which es- 
sentially reflects the late formation of high-mass clusters 
in ACDM (e.g., Sadeh & Rephaeli 2008). At z = the 
mean beta profile of high-mass halos appears to be some- 
what steeper at r > 0.4r V j r than that of low-mass halos, 
in agreement with the fact that for the same radii it 
tends to higher levels at high redshifts, as illustrated in 
the upper panel of figure |3l 

As mentioned in § [4J the criteria according to which 
halos are classified as relaxed or unrelaxed are quite ar- 
bitrary. We chose in our analysis to compare among the j3 
profiles by setting two comparable numbers correspond- 
ing to the most and least relaxed phases. In figure 2] we 
plot the velocity anisotropy profile of relaxed versus unre- 
laxed halos using spherical shells when the distinction is 
made according to the s a f f se t criterion. Using the virial 
ratio to distinguish between relaxed and unrelaxed ha- 
los gave a very similar ft profiles, and therefore these are 
not shown here. At radii smaller than the virial radius 
applying the s ffset criterion results in flattened velocity 
anisotropy profiles of the unrelaxed halos with respect to 
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Fig. 3. — Velocity anisotropy profiles of halos with different 
masses. Top panel: profiles of the 100 most (solid curve) and least 
(dashed curve) massive halos at z = 0. Bottom panel: profiles of 
the 10 most (solid) and least (dashed) massive halos at z = 2. 


















s 






-1 




\ 

\ 


\ V /\ 


-2 
-3 






\'\ \ 






Relxed: s „ <0.01 5, 267 halos 
offset 


> 1 


-4 




Unrelxed: s^. >0.17, 264 halos 


> 1 
1 1 
'l 








V 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Radius [r .1 



Fig. 4. — Velocity anisotropy profiles of relaxed (blue solid curves) 
and unrelaxed (red dashed curves) halos analyzed in spherical 
shells. Relaxation gauged by the s ff aet criterion: relaxed and 
unrelaxed halos have s fr se t < 0.015 (264 halos) and s ffsot > 0.17 
(267 halos), respectively. As in the previous plots, the central curve 
represents the mean value, whereas the upper and lower curves de- 
scribe the ± 1— <7 uncertainty range. 
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Fig. 5. — Velocity anisotropy profiles when the shells are de- 
scribed as spherical (red dash curves) or elliptical (blue solid 
curves). Note, the radius in the elliptical case is the semi- major 
axis. 



Fig. 7. — The fraction of relaxed halos at different ellipticities 
and for different relaxation criteria, 2T/\U\ = 1.5 (blue circles) 
and s ff se t = 0.05 (red asterisks). The uncertainties were taken 
to be Poisson statistics, i.e. AJV = 1/yN. 
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Fig. 6. — Velocity anisotropy profiles in elliptical shells. The 
average ft profile of spherical halos, q > 0.95, (blue solid curves) is 
compared to a elliptical halos, 0.5 < q < 0.517, (red dash curves). 

the relaxed halos. 

To assess the impact of aspehrical halo configuration, 
we plot in figure [5] the velocity anisotropy profiles in 
both spherical and elliptical shells for all high-mass ha- 
los (M v i r ^ 10 14 hp \ Mq). It is interesting to note that 
the average (3 profile is almost the same in spherical and 
elliptical shells till ~ (0.7 — l)r v ; r . From these region, 
however, the scatter of /3 using elliptical shells is much 
smaller. 

In figure [6] we used elliptical shells. We derived the 
average (3 profile of spherical halos, q > 0.95, (blue solid 
curves) and compared it to the one derived from elliptical 
halos, 0.5 < q < 0.517, (red dash curves). The j3 values 
of spherical halos are lower than the ones of elliptical 
halos. In addition, the scatter in /3 is larger in spherical 
halos. 

In figure [7] we plotted the connection between halo 
cllipticity and relaxedness. We took only halos with 
q > 0.4 since more elliptical halos are very few and 
therefore give poor statistics. In addition, for halos with 
q < 0.4 the values of both relaxation criteria are strongly 
dependent on the length of the ellipse major axis, which 
is likely due to the fact that many of them are in the 
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Fig. 8. — Velocity anisotropy profiles of relaxed (blue solid curves) 
and unrelaxed (red dashed curves) halos in elliptical shells. Relax- 
ation gauged by the s a f j ae t criterion: relaxed and unrelaxed halos 
have s g set < 0.013 (311 halos) and s ff set > 0.13 (323 halos), re- 
spectively. As in the previous plots, the central curve represents 
the mean value, whereas the upper and lower curves describe the 
± 1— a uncertainty range. 

process of a major merger and highly unrelaxed. The 
threshold was chosen so the two will have about the same 
normalization 

The (3 profile of relaxed versus unrelaxed halos in el- 
liptical shells are not appreciably different than the ones 
using spherical shells. In figure [8] we plot the velocity 
anisotropy profile of relaxed versus unrelaxed halos using 
elliptical shells when the distinction is made according to 
the s f / se t criterion. The profiles are similar to those for 
spherical shells except for a smaller decline and with a 
smaller scatter at large radii. 

5.2. 7-/? ratio 

As was mentioned in § [TJ the question of whether 7 
and f3 are correlated is of both theoretical and practi- 
cal interest. In figure |H] we show the velocity anisotropy 
vs. the radial density slope for all shells in all halos (left 
panel), all shells of relaxed halos according to the virial 
relation criterion 2T/\U\ < 1.35 (middle panel), and all 



G 



shells of highly relaxed halos with 2T/\U\ < 1.35 and 
^offset < 0.025 (right panel). For each halo we checked the 
maximum grid level and determined the halo minimum 
spatial resolution. Unresolved shells were not included in 
the analysis. The black curve reproduces the Hansen & 
Moore relations for the —4 < 7 < range. In figurefTUlwc 
drew the same quantities for all shells in all halos, assum- 
ing NPW-distributed density profiles. Note that in this 
plot the 7 range is (~ —2.9, —1), not (—3,-1), due to the 
finite binned values of the radius. Finally, figure [TT] de- 
scribes the velocity anisotropy against the radial density 
slope of the four inner (0 < r < 0.3r v j r , top panel) and all 
the other (0.3r v i r < r < r v i T , bottom panel) shells. This 
boundary value was chosen since shells included within 
this radius, 0.3r v ; r , display the strongest 7-/? correlation. 

We repeated this above analysis with triaxial halos; 
the results were essentially the same as those obtained 
for spherical halos. 

6. SUMMARY 

Significant progress has recently been made in the abil- 
ity to deduce the kinematic properties of DM in clusters 
from galactic dynamics and X-ray measurements. Com- 
parison of these properties with results from numerical 
simulations can clearly test analysis methods and add 
new insights on DM phase space occupation. We pre- 
sented results from an analysis of DM velocities in 6019 
halos with masses M ^ 10 14 hg j. M at redshift 2 = 0, 
drawn from one of the largest ever hydrodynamic cosmo- 
logical AMR simulations. 

Our study indicates that the profiles of cluster DM 
velocity dispersions have a similar pattern for all halo 
masses, redshifts, and relaxation stages, even though 
there is a considerable scatter in magnitude due to the 
large differences in the j3 profile of individual halos. A 
typical behavior is a rising /3 profile from a nearly vanish- 
ing central value, leveling off at r ~ 0.2r V j r , out to large 
radii of at least (1.5 — 2)r v j r . Lower mass halos at lower 
redshifts have, on average, lower j3 levels at r = r v ; r , and 
therefore even shallower profiles. For example, low mass 
halos of M f=s 10 14 hp \ M Q at z — attain levels of 
j3(r — r V i r ) « 0.2, while high mass halos, i.e. M > 10 15 
h^y M , reach slightly higher values, f}(r — r v i r ) w 0.4. 
This behavior could possibly be due to the ability to 
reach higher accretion velocities at lower redshifts when 
the background density is much lower. Lau, Nagai, & 
Kravtsov (2010) showed that the inclusion of radiative 
cooling and star formation in the simulation slightly low- 
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ers the j3 values at z = 0. At higher redshifts, e.g. z — 2, 
we found no significant difference between the (3 profile of 
high and low mass halos, even though the scatter is lower 
than at z = 0. This is also in agreement with the trend 
found by Lau, Nagai, & Kravtsov (2010). When gauged 
by s f f se t , very relaxed halos have practically similar /3 
levels to the corresponding levels of low mass halos, i.e. 
(3{r = r v j r ) s=a 0.2 — 0.3 at z = 0. However, when gauged 
by the virial ratio, the very relaxed and unrelaxed halos 
have similar /3 profile. 

The shapes of cluster halos are of considerable inter- 
est that has recently been enhanced by observational ca- 
pabilities (e.g. the CLASH project). Observed samples 
may, however, not be large enough for a conclusive exam- 
ination. We show here that cluster shape measurements 
can be combined together with j3 values at large radii 
for stronger constraints. It is important to measure (3 at 
large radii, > 0.3r V j r , since on the one hand it can be de- 
rived with a high level of confidence (Lemze et al. 2011); 
on the other hand, at these radii the j3 values tend to be 
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halos. Note that this difference depends on the value set 
for the major axis. It is also important to analyze (3 pro- 
files in elliptical shells as this significantly reduces their 
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We have employed two different relaxation criteria and 
found that more spherical halos tend to be more relaxed. 
This can be explained by the evolution of clusters from 
highly aligned and elongated systems at early times to 
lower alignment and elongation at present, which reflects 
the hierarchical and filamentary nature of structure for- 
mation (Hopkins, Bahcall, & Bode 2005). Indeed, Vera- 
Ciro et al. (2011), who analyzed Aquarius data, found 
that q increases with time. 

Lastly, we find that there is some correlation between 
7 and f3 at low radii, r < 0.3r v j r , and that such a cor- 
relation can be induced at all radii merely by assuming 
a prescribed DM density profile. The level of 7 - j3 cor- 
relation is very low at large radii, r > 0.3r v ; r , even for 
very relaxed halos. Repeating the same analysis with 
elliptical shells led to the same result. 
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Fig. 10. — Velocity anisotropy vs. radial density slope at z = 0. 
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NFW fit (blue dots). The HM06 relation is plotted in the -4 < 
7 < range (black solid curve). 
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Fig. 11. — Velocity anisotropy vs. radial density slope at z = 0. 
Top panel: plotted are values for the four inner shells, < r < 
0.3r v ; r . Bottom panel: values for the outer shells, 0.3r v ; r < r < 
r v ; r . The HM06 relation is depicted in the —4 < 7 < range (black 
solid curve). The vertical line at 7 ~ —2 is an artifact due to the 
fact that 7 here is a discrete slope profile calculated between two 
bins with a low number of particles. This artifact can also be seen 
in fig. [9] though it is more prominent here. 



